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combined to substantiate our claims. 
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1. Introduction 

Sociophysics is a long standing^ research field addressing issues related to 
the characterization of the collective social behavior of individuals, such 
as culture dissemination, the spreading of linguistic conventions, and the 
dynamics of opinion formation. ^"^ These are all interdisciplinary applica- 
tions which sit the interface of different domains. The challenge is in fact to 
model the dynamical evolution of an ensemble made of interacting, micro- 
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constituents and infer the emergence of collective, macroscopic behaviors 
that are then eventually accessible for direct experimental inspection. Agent 
based computational models are widely employed in sociophysics applica- 
tions and also adopted in this paper. They provide in fact a suitable setting 
to define local rules which govern the evolution of the microscopic con- 
stituents. 

In recent years, much effort has been devoted to the investigation of social 
networks, emerging from interaction among humans. In the sociological lit- 
erature a main distinction has been drawn between smalF and large*^"'^^ 
groups, as depending on its intrinsic size, the system apparently exhibits 
distinct social behaviors. Up to a number of participants of the order of a 
dozen, a group is considered small. All members have a clear perception 
of other participants, regarded as individual entities: Information hence 
flows because of mutual relationships. Above this reference threshold, se- 
lected individuals see the vast majority of the group as part of a uniform 
mass: There is no perception of the individual personality, and only aver- 
age behaviors matter. This distinction is motivated by the fact that usually 
humans act following certain prebuilt "schemes" ^^'^^ resulting from past 
experiences, which enables for rapid decision making whitout having to 
necessarily screen all available data. This innate data analysis process al- 
lows one for a dramatic saving of cognitive resources. 

These conclusions have been reached on the basis of empirical and qualita- 
tive evidences. ^'^"^"^^ We are here interested in detecting the emergence of 
similar phenomena in a simple model of opinion dynamics. ^®'^ As we shall 
see in our proposed formulation agents posses a continuous opinion on a 
given subject and possibly modify their beliefs as a consequence of binary 
encounters. 

The paper is organized as follows. In the next section the model is intro- 
duced. Forthcoming sections are devoted to characterizing the quantities 
being inspected and develop the mathematical treatment. In the final sec- 
tion we sum up and comment about future perspectives. 

2. The model 

We shall investigate the aforementioned effects related to the size of the 
group of interacting individuals (social group), within a simple opinion dy- 
namics model, recently introduced in^® . For the sake of completeness we 
hereafter recall the main ingredients characterizing the model. The inter- 
ested reader can refer to^^ for a more detailed account. 
We consider a closed group composed by N individuals, whose opinion on a 
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given issue is represented by a continuous variable Oi , scanning the interval 
[0, 1]; moreover each agent i is also characterized by the so-called affinity, 
a real valued vector of iV — 1 elements, labeled aij, which measures the 
quality of the relationship between i and any other actor j belonging to the 
community. 

Agents interact via binary encounters possibly updating their opinion and 
relative affinity, which thus evolve in time. Once agents i and j interact, via 
the mechanism described below, they converge to the mean opinion value, 
provided their mutual affinity scores falls below a pre-assigned threshold 
quantified via the parameter Oc- In formulae: 

Ol+' ^Ol-\ I\0% Ti (4) & Of 1 ^ O*- - ^ ^^% Ti (a*,) , (1) 

where AO*^ = 0\ - O) and Ti [x) = \ [tanh(/3i(a; - a^)) + 1]. The latter 
works as an activating junction defining the region of trust for effective 
social interactions. On the other hand bearing close enough opinions on 
a selected topic, might induce an enhancement of the relative affinity, an 
effect which is here modeled as: 

a*+i=a^+4(l-4)r2(AOy & a*f = a*,+4(l-4)r2(AO*,) 

(2) 
being r2 (x) = — tanh(/32(|a;| — AOc))- This additional activating function 
quantifies in AOc the largest difference in opinion (AO*) which yields to 
a positive increase of the affinity amount a* . The parameters /3i and P2 
are chosen large enough so that Fi and r2 are virtually behaving as step 
functions. Within this working assumption, the function Fi assumes value 
or 1, while F2 is alternatively — f or +1, depending on the value of their 
arguments ^. 

The affinity variable, a*,, schematically accounts for a large number of hid- 
den traits (diversity,personality, attitudes, beliefs...), which are nevertheless 
non trivially integrated as an abstract simplified form into the model. Note 
also that the affinity accounts for a memory mechanism: indeed once two 
agents meet, the outcome of the interaction in part depends on their history 
via the affinity scores. 
To complete the description of the model let us review the selection rule 



''We shall also emphasize that the logistic contribution entering Eq. (2) maximizes the 
change in affinity when a*. Ki 0.5, corresponding to agents i which have not yet build 
a definite judgment on the selected interlocutor j. Conversely, when the affinity is close 
to the boundaries of the allowed domain, marking a clear view on the worth of the 
interlocutor, the value of a* ■ is more resistant to subsequent adjustments. 
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here implemented. Each time step a first agent i, is randomly extracted, 
with uniform probability. Then a second agent j is selected, which mini- 
mizes the social metric D\, and time t. This is a quantity defined as: 

D%^d%+M,{Q,a), (3) 

where d* — |A0* |(1 — a* ) is the so-called social distance and Mj{Q,cf) 
represents a normally distributed noise with zero mean and variance a, 
that can be termed social temperature}^ The rationale inspiring the mech- 
anisms here postulated goes as follows: The natural tendency for agent i to 
pair with her/his closest homologous belonging to the community (higher 
affinity, smaller opinion distance), is perturbed by a stochastic disturbance, 
which is intrinsic to the social environment (degree of mixing of the popu- 
lation) . 

The model exhibits an highly non linear dependence on the involved pa- 
rameters, oic, AOc and <J. In a previous work^® the asymptotic behavior of 
the opinions dynamics was studied and the existence of a phase transition 
between a consensus state and a polarized one demonstrated. It should 
be remarked however that the fragmented case might be metastable; in 
fact if the mean separation between the adjacent opinion peaks is smaller 
than the opinion interaction threshold, AOc, there always exists a finite, 
though small, probability of selecting two individuals belonging to different 
clusters, hence producing a gradual increase in the mutual affinities, which 
eventually lead to a merging of the, previously, separated clusters. This fi- 
nal state will be achieved on extremely long time scales, diverging with the 
group size: socially relevant dynamics are hence likely to correspond to the 
metastable regimes. 

A typical run for A^ = 100 agents is reported in the main panel of Fig. 1, 
for a choice of the parameters which yields to a monoclustercd phase. This 
is the setting that we shall be focusing on in the forthcoming discussion: 
Initial opinions are uniformly distributed in the interval [0, 1], while a^- are 
randomly assigned in [0, 1/2] with uniform distribution, parameters have 
been fixed to ac = 0.5, AOc = 0.5 and a = 0.01. 

Once the cluster is formed, one can define the opinion convergence time, 
Tc, i.e. time needed to aggregate all the agents to the main opinion group ^ . 
A second quantity Tq can be introduced, which measures the time scale for 
the convergence of the mean group affinity to its asymptotic value 1. The 
latter will be rigorously established in the next section. 



''We assume that a group is formed, i.e. aggregated, once the largest difference between 
opinions of agents inside the group is smaller than a threshold, here 10~*. 
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Such quantities are monitored as function of time and results are schemati- 
cally reported in Fig. 1. As clearly depicted, the evolution occurs on sensibly 
different time scales, the opinion converging much faster for the set of pa- 
rameters here employed. 

In the remaining part of this paper, we will be concerned with analyzing the 
detail of this phenomenon highlighting the existence of different regimes as 
function of the amount of simulated individuals. More specifically, we shall 
argue that in small groups, the mean affinity converges faster than opinions, 
while the opposite holds in a large community setting. Our findings are to 
be benchmarked with the empirical evidences, as reported in the psycholog- 
ical literature. It is in fact widely recognized that the dynamics of a small 
group (workgroup) proceeds in a two stages fashion: First one learns about 
colleagues to evaluate their trustability, and only subsequently weight their 
input to form the basis for decision making. At variance, in large communi- 
ties, only a few binary interactions are possible among selected participants 
within a reasonable time span. It is hence highly inefficient to wait accu- 
mulating the large number of information that would eventually enable to 
assess the reliability of the interlocutors. The optimal strategy in this lat- 
ter case is necessarily (and solely) bound to estimating the difference in 
opinion, on the topic being debated. 




Fig. 1. Opinions as function of time. The underlying network is displayed at different 
times, testifying on the natural tendency to evolve towards a coherent ensemble of afBne 
individuals. Tc and Ta are measured according to our conventions. 
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3. The social network of affinities 

In our model the affinity enters the selection mechanism that makes agents 
to interact. We can restate this fact by postulating the existence of an un- 
derlying social network, which drives the interactions and thus the opinion 
flow. In this perspective, the affinity can be seen as the adjacency matrix of 
a weighted '^ graph. In such a way we arc formally dealing with an adaptive 
social network^*'^^ : The network topology influences the opinion dynamics, 
the latter providing a feedback on the network itself. In other words, the 
evolution of the topology is inherent to the dynamics of the model because 
of the proposed self-consistent formulation and not imposed a priori as an 
additional, external ingredient, i.e. rewire and/or add/remove links accord- 
ing to a given probability^°'^^ once the state variables have been updated. 
From this point of view, the mean group affinity is the averaged outgoing 
degree - called for short "the degree" in the following - of the network: 

<^>w = ^E^*' (4) 

i 

where the degree of the i-th node is fc* = ^ a\J{N — 1). The normalizing 
factor A^ — 1 allows for a direct comparison of networks made of a different 
number of agents. Let us observe that we chose to normalize with respect 
to A^ — 1 because no self-interaction is allowed for. 

In the left panel of Fig. 2 we report the probability distribution function for 
the degree, as a function of time. Let us observe that the initial distribu- 
tion is correctly centered around the value 1/4, due to the specificity of the 
chosen initial condition. The approximate Gaussian shape results from a 
straightforward application of the Central Limit Theorem to the variables 
(fc*)i=i^...^Ar. In the right panel of Fig. 2 the time evolution of the mean 
degree < k > (t) is reported. Starting from the initial value 1/4, the mean 
degree increases and eventually reaches the value 1, characteristic of a com- 
plete graph. As previously mentioned this corresponds to a social network 
where agents are (highly) affine to each other. In the same panel we also 
plot the analytical curve for < /c > (t), as it is determined hereafter. 

From the previous observation, it is clear that the time of equilibration 
oi < k > (t) provides an indirect measure of the convergence time for 
a* ■ -^ 1. The affinity convergence time, Tq can be thus defined as: 

T„ = min{i >0:< k> (t)>'q}, (5) 



'^In fact the trustability relation is measured in terms of the "weights" a* G [0, 1] 
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Fig. 2. Time evolution of the degree probability distribution function and the averaged 
degree. Left panel : Several histograms representing successive snapshots of the dynamics 
are displayed: t = (black online), two generic intermediate times (wight), t = Tc (red 
online) and t = Ta (blue online) ; histograms are normalized to unity. Right panel : < k > 
versus time. Symbols refer to direct simulations. The solid lines are obtained as a one 
parameter fit of the theoretical expression Eq. (10). 

where rj € (0, 1) is a threshold quantity. The closer 77 is to 1 the stronger the 
condition on the value of a* (the larger the value for Ta) which identifies 
an affine unit. In the following we shall assume rj = 3/4. 



4. Results 

The aim of this section is to analyze the behavior of Ta and of < A: > (t) 
so to provide an analytical interpretation for the results presented in the 
previous section. 

From the definition of mean degree Eq. (4) one can compute the time 
evolution of < fc > (i) in a continuous time setting: 



d < k > 
Jt 



1 



da* 



N{N - 



1)^ 



dt 



(6) 



To deal with the update rule for the affinity Eq. (2) we assume that we can 
decouple the opinions and affinity dynamics by formally replacing the acti- 
vating function r2(AO*.) with a suitably tuned constant. Strictly speaking, 
this assumption is correct only for AOc = 1, in which case the opinions are 
always close enough so to magnify the mutual affinity scores of the interact- 
ing agents, as a results of a self consistent evolution. For the general case, 
AOc < 1, what we are here proposing is to replace r2(A0,-) by some ef- 
fective value termed jeff , determined by the dynamics itself. The technical 
development that yields to the optimal estimate of "f^ff will be dealt in the 
appendix. 
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Under this assumption, the equation Eq. (2) for the evolution of a* admits 
the following continuous version : 

-^ = 7e//a*, (1 - ay , (7) 

which combined to equation Eq. (6) returns the following equation for the 
mean degree evolution : 

d<k > 



dt 



= 7ei/(<fc>(i)-<(4)2>) . (8) 



Assuming the standard deviation of (a*,) to be small "^ for all t, implies 
< {a\j Y >^< Oiij >^~< k >^, which allows us to cast the previous 
equation for < fc > in the closed form : 

^^=7e//(<fc>-<fc>^). (9) 

This equation can be straightforwardly solved to give: 

ko + {l- ko)e ^'f'^ 

where fco =< k > (0). We can observe that such solution provides the 
correct asymptotic value for large t. Moreover j^ff plays the role of a char- 
acteristic time and in turn enables us to quantify the convergence time Tq 
of the affinity via : 

T^.^U^il-M).^. (U) 

leff \ko{l-T])J Jeff 

In the appendix we determine ° the following relation which allows to ex- 
press Jeff as a function of the relevant variables and parameters of the 
models, i.e. Tg, AOc and N: 

^^ff = ^ + ^P^ (12) 

where p = -(1 + 2AOclog(AOc) - ^Oj), thus recalling Eq. (11) we can 
finally get: 

^" - ^'^' (l - W) ■ (13) 



'^This assumption is supported by numerical simulations not reported here and by the 

analytical considerations presented in.'^'^ 

°This is case a) of Eq. (A. 7). In the second case the result is straightforward: Tq = 
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From previous works^^'^^ we know that the dependence of Tc on N , for 
large N , is less than quadratic. Hence, for large N , the second term in 
the parenthesis drops, and we hence conclude that, the affinity convergence 
time grows like T^ ^ N"^, as clearly shown in the main panel of Fig. 3. 
The prefactor'e estimate is also approximately correct, as discussed in the 
caption of Fig. 3. 

The above results inspire a series of intriguing observation. First, it is im- 
plied that the larger the group size the bigger Ta with respect to Tc- On 
the contrary, making N smaller the gap progressively fades off. Dedicated 
numerical simulations (sec left inset of Fig. 3) allows to indentify a turning 
point which is reached for small enough values of N: there the behavior 
is indeed opposite, and, interestingly, Tc > Ta. The transition here repro- 
duced could relate to the intrinsic peculiarities of the so called "small group 
dynamics" to which we made reference in the introductory section''^^*' . Fur- 
thermore, it should be stressed that the critical group size determining the 
switching between the two regimes here identified, is quantified in A'^ ~ 20, 
a value which is surprisingly closed to the one being quoted in social studies. 



• 
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= 0.4 
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= 0.6 



1 1 4 * ^ 




20 40 60 SO 100 



."iO 100 150 200 liO .100 350 400 450 500 
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Fig. 3. Main panel: Ta/N^ vs N for different values of the parameter AOc. The data 
approach a constant value {Ta/N^ ~ 1.72) clearly indicating that the time of convergence 
of the affinity matrix scales quadratically with the number of agents, in agreement with 
the theory. The asymptotic value estimated by our theory is 2.19, the discrepancy being 
therefore quantified in about 15 %. Left inset: T^/N'^ and T^/N'^ vs N for AOc = 0.5. 
As predicted by the theory and the numerics a crossover is found for groups for which 
opinions converge slower than the affinities: this is the signature of a distinctive difference 
in the behavior of small and large groups, numerically we found that this difference is 
effective for A'^ ~ 20. Right inset: < fc > vs t/N^ is plotted for two different values of N . 
As expected the two curves nicely collapse together. 
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5. Conclusion 

In this paper we study a model of continuous opinions dynamics already 
proposed in^^ , which incorporates as main ingredient the affinity between 
agents both acting on the selection rule for the binary interactions as well 
entering the postulated mechanism for the update of the individual opin- 
ions. 

Analyzing the model in the framework of adaptive networks we have been 
able to show that the sociological distinction between large and small groups 
can be seen as dynamical effect which spontaneously arises in our system. 
We have in fact proven that for a set of realistic parameters there exists 
a critical group size, which is surprisingly similar to the one reported in 
the psychological literature. Below this critical value agents first converge 
in mutual affinity and only subsequently achieved a final consensus on the 
debated issue. At variance, in large groups the opposite holds : The con- 
vergence in opinion is the driving force for the aggregation process, affinity 
converging on larger time scales. 

Appendix A. Computation of 7e// 

The aim of this paragraph is to provide the necessary steps to decouple 
the opinion and affinity dynamics, by computing an effective value of the 
activating function r2(AO*^), hereby called "feff- This will be obtained by 
first averaging r2 with respect to the opinions and then taking the time 
average of the resulting function: 

1 /■* 

7e// = lim - / dT <T>op (t) , (A.l) 

where the opinion-average is defined by: 

<r>„p(t)= [ dx [ dyT2{\x-y\)f{x)f{y), 
Jo Jo 

being /(•) the opinions probability distribution function. To perform the 

computation we will assume that for each t, the opinions are uniformly 

distributed in the interval [a{t), a{t) + L{t)] where L{t) ~ 1 — t/Tc and Tc 

is the opinion convergence time, hence /(•) = l/(NL{t)). This assumption 

is motivated by the "triangle-like" convergence pattern as clearly dsplayed 

in the main panel of Fig. 1. 

Assuming (32 large enough, we can replace r2 by a step function. Hence: 

2 fa{t)+L(t) /.a(t)+L{t) 



<T>op{t) ^ ^,.,,^,,, / dx I dyx{x,y), (A.2) 
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x{x,y) 



where x(a;, y) is defined by (see also Fig. Al) 

' 1 if |a; - y| > AOc, i.e. in tlie triangles Ti U Ta = Q \ Z? 
— 1 otherwise, i.e. in D , 

(A.3) 
where Q is the square [a, a + L] x [a,a + L]. 

Let us observe that this applies only when L{t) > AOc (see left panel of 
Fig. Al); while if L(t) < AOc the whole integration domain, [a, a + L] x 
[a,a + L], is contained into the ja; — j/| < AOc (see right panel of Fig. Al). In 
this latter case, the integration turns out to be simpler. In other words the 



y=x4-AO^ 



a(t)+L(tl 




Fig. Al. The geometry of the integration domains. On the left panel the case L > AOc 
while on the right one the case L < AOc- 



integral in Eq. (A. 2) corresponds to measure the area, \D\, of the domains 

shown in Fig. Al with a sign. Let us perform this computation according 

to the two cases: L > AOc or L < AOc- 

In the case L > AOc, the area of D is given by \D\ = \Q\ - |Ti| - IT2I, 

hence 

1 



< r >op (t) = - 



1 



1 



H^l + I^il 



T2 = - 



-L' 



, {L - AOc) 



1 



(-101 



1-2 1- 



2\Ti\+2\T2 
AO. 



L 



On the other hand if L < AOc, because the square Q is completely con- 
tained into the domain |a; — y | < AOc where % is equal to —1, we easily get: 
< r >op (t) = —jfTjj{—L^) = 7^, if i < AOc- This last relation together 
with Eq. (A. 4), can be castcd in a single formula: 



<r> 



1 



1-2 1- 



AOr 



e (L - AOc) 



(A.5) 



(A.4) 
(if L > AOc) ■ 
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where is the Heaviside function, <d{x) — 1 ii x > and zero otherwise. 
To conclude we need to compute the time average of < F >op {t). Using 
once again the "triangle-hke" convergence assumption for the opinions, i.e. 
L(t) = 1 — t/Tc, where Tc is the opinion convergence time, we get: 



7e// 



1 


/•* 1 


hm — 


/ dr , 


t^oo t 


N^ 






(A.6) 



This integral can be explicitly solved to give: 



ifr„<T,, 



p+i 

Ar2 



2^ 



where p=-{l + 2AOclog(AOc) - AO, 
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